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We investigate the robustness of the microscopic reversibility in open quantum systems which is 
discussed by Monnai [arXiv:1 106.1982 (2011)]. We derive an exact relation between the forward 
transition probability and the reversed transition probability in the case of a general measurement 
basis. We show that the microscopic reversibility acquires some corrections in general and discuss 
the physical meaning of the corrections. Under certain processes, some of the correction terms 
vanish and we numerically confirmed that the remaining correction term becomes negligible; for 
such processes, the microscopic reversibility almost holds even when the local system cannot be 
regarded as macroscopic. 



I. INTRODUCTION 



Understanding of non-equilibrium quantum dynamics has been eagerly pursued since many experiments are done 
under non-equilibrium situations and should be also treated quantum mechanically. While there are many physical 
quantities and relations to characterize the properties of the equilibrium systems, very few are known for the systems 
out of equilibrium. The linear response theory, which describes the non-equilibrium quantities in terms of the equilib- 
rium quantities, is a powerful tool to investigate such systems. It is, however, restricted to the case where the systems 
are close to equilibrium. In order to describe the behaviors of strongly non-equilibrium systems, it may be required 
to find some relations which contain detailed information of the dynamics. 

For exploring the character of the strongly non-equilibrium dynamics of quantum systems, symmetry relations for 
non-equilibrium states such as fluctuation theorems are of great importance [lH9|. In particular, significant attentions 
are paid to the relation between the transition probabilities of forward processes and the transition probabilities of the 
corresponding reversed processes 0, H, H-El . Such relations are appealing particularly because they are not restricted 
to the close-to-equilibrium states. As the forward process of the fluctuation theorems, it is typical to consider a 
driven system which is in the thermal equilibrium state at t = and then is controlled by a time-dependent external 
parameter X(t). We consider the case where the states at the beginning and the ending of the process are measured, 
and thus the states are projected. We can interpret that the state hopped from the one to the other; the transition 
probability of the forward process is defined as the probability that such a transition occurs. For the reversed process, 
it is common to consider the following in quantum systems jH S HI, [H, Hj] . The initial density matrix of the reversed 
process is set to be the same as that of the forward process but with the value of the external parameter A(T), where 
T is the final moment of the forward process. Then, the system is driven with a time-reversed protocol A(T — i); 
the transition probability of the reversed process is defined as the probability that the opposite transition occurs 
compared to the transition of the forward process. 

They are the typical settings for the study of fluctuation theorems. From the viewpoint of quantum operation, 
however, it may be natural to set the time-evolved (generally non-equilibrium) state as the initial density matrix of 
the reversed process. 

As another point for the study of open quantum systems, the choice of the measurement process is especially 
significant; measurement processes on a reservoir Q or both on a local system and a reservoir @, H, Q are often 
considered. Nevertheless, sometimes it is more natural to consider the measurement solely on the local system; in the 
present paper, we consider such a case. There are fluctuation theorems of open quantum systems which are written 
solely in the terms of local system though they do not take into account of the measurement process in order to 
discuss the forward and the reversed processes [T^, EH . 

The above two points were taken into account in a recent study by Monnai [l6jj. The study pointed out the 
significance of the microscopic reversibility in open quantum systems as a kind of symmetry relation similar to but 
different from the fluctuation theorems [16[ ■ It defines the reversed process as the one from the time-evolved state 
and considers the measurements on the local system only. 

The discussion in Ref. [l(| for open quantum systems is as follows. The total Hamiltonian consists of the Hamiltonian 
of a local system H s (X(t)) which is controled by external forces with the parameter X(t), the Hamiltonian of a reservoir 
H r , and the Hamiltonian of coupling between them H c , i.e. , 



fltot(t) = 4(A(t)) + H T + H c . 



(1) 
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Let us consider the process where we measure the states of the local system at t = and t = T. Thoughout this paper, 
we employ the Schrodingcr picture and only consider the projection measurement as the measurement protocol. The 
measurement basis at t — can be different from the one at t = T. We refer to the measured states as \n(Q)) and 
\m(T)) and to the probability of such a transition as p-p{\n{Q)) — > \m{T))). Note that the time variables of |n(0)) and 
|m(T)) merely indicate the moments that the measurements are done along the forward process; they do not mean 
that those measurement bases are time dependent. Next, the reversed process is defined as follows; as the initial 
state, we prepare the state p(T) that evolved from t = to t = T without the measurement at t = 0, and then drive 
the system from t = T to t = IT with the time- reversed protocol A(2T — t). The reversed transition probability is 
defined to be the probability of observing Q\m(T)) at t = T and 0|n(O)) at t — 2T under the time-reversed process, 
where is the time-reversal operator. We refer to such a transition probability as pR,(0|m(T)) —> 0|n(O))). Monnai 
then showed the equality [l6[ 

pp(|n(0)) -> |m(T))) = p R (e\m(T)) -> 9|n(0)» (2) 

under the following conditions: (i) the total system is a product state at t = 0; (ii) the local system is macroscopic, 
so that the contribution from the coupling Hamiltonian H c is extremely small compared to the ones from the local 
system H s and the reservoir J? r , and we thereby have U {T){yJ p s (0) <£> \J /5 r (0))C/ t (T) ~ y/ p s (T) ® y/p T (0), where p s 
and p r are the density matrices of the local system and the reservoir, respectively; and (iii) the measured states at 
t = and t = T are the eigenstates of the density matrix of the local system. We call the equality ^ the microscopic 
reversibility in open quantum systems. Note that it is a relation about the local system; no measurements are done 
on the reservoir. 

The main purpose of this paper is to investigate the robustness of the microscopic reversibility ([2"]). Although we 
assume that the initial state is a product state, we will allow the final state to be arbitrary; we will not assume 
the local system to be macroscopic and we will consider an arbitrary measurement bases. We will also assume that 

(m(T)\Q p{T)e\m(T)) = (m{T)\p(T)\m{T)); the probability to obtain the resulting state |m(T)) is equal to that of 
the time-reversed state. As a result, we will show that the microscopic reversibility does not hold exactly in general; 
it acquires correction terms. The origin of the corrections is the effect that the measurement processes destroy the 
quantum coherence of the system that we measure. Although the microscopic reversibility is broken in general, if we 
measure the eigenstatc of the density matrix of the local system at t = 0, the form of the correction becomes very 
simple. In the case of a thermal relaxation process, we numerically confirmed that the correction term is small enough 
compared to the forward and the reversed transition probabilities in open quantum systems. 

This paper is organized as follows: In Sec. [TTl we will derive the microscopic reversibility in isolated quantum 
systems with correction terms under general measurement bases. As a simple example, we will consider the case of a 
free particle; we will show that a correction term can be very large in this case, and thus we cannot see the microscopic 
reversibility at all. In Sec. Mil wc will derive the microscopic reversibility of open systems with correction terms. If 
the initial and the final states are product states and the eigenstates of the density matrices of the local system are 
measured, we will show that the correction terms vanish and the microscopic reversibility holds exactly. We will also 
show that, if we disconnect the local system from the reservoir, the correction term then becomes constant. In Sec. lIVl 
we analyze the details of the corrections. Finally, in Sec. [Vj we numerically compute a correction to the microscopic 
reversibility for a one-dimensional spin chain; we regard the first two spins as the local system and the rest as the 
reservoir. The result shows that the correction is relatively small, so that the microscopic reversibility almost holds 
even when the local system cannot be regarded as macroscopic. 

II. MICROSCOPIC REVERSIBILITY IN ISOLATED SYSTEMS 

We first describe the microscopic reversibility in isolated quantum systems with the same notations and processes 
as in Sec. |U The forward and the reversed transition probabilities read [l| 

MHO)) -> |m(T)» = (m(T)|^|n(0)>(n(0)|p(0)|n(0)>(n(0)|^|m(T)>, (3) 
p R (e\m(T)) -> e|n(0))) = (n(0)|^t| m (T))(m(T)|ep(T)6|m(T))(m(T)|^|n(0)). (4) 

<r 

Throughout this paper, we consider the case where (m(T)\Q p(T)Q\m(T)) = (m(T)\p(T)\m(T)}. This condition is 
satisfied, for example, in the case where the states or the total Hamiltonian are invariant under the time reversal. In 
Sec. Ill Al we treat the former case, the measurement of the position of a particle. In Sec. [V] we treat the latter case, 
the spin system which is invariant under the flip of all the spins. Then we have the following relation between the 
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forward transition probability and the reversed transition probability: 

MHO)) -> |m(T))) = f nm (T,0) PR (e\m(T)) 



6|n(0))), 



(5) 



where / nm (T,0) = (n(0)|p(0)|n(0)}/(m(T)|p(T)|m(T)}. Let us consider the trivial case where we choose the measure- 
ment bases at t = and t = T as the eigenstates of the density matrices p{0) and p(T), respectively. If we denote 
|n'(T)) := £/|n(0)), we have (m(T)\U\n(0)) = 5 n > m , and therefore 



(n(0)|#0)|n(0)} = (n(0)\U^Um^UH0)) = Sn',n(m(T)\p(T)\m(T)). 
Hence the microscopic reversibility trivially holds: 

MHO)) -> |m(T))) = p a (6|m(T)) -> 9|n(0)». 



(6) 



(7) 



In many cases, however, it is difficult to detect the eigenstates of an arbitrary density matrix or make the system 
have a density matrix whose eigenstates coincide with the measurement basis that we choose. In the case where the 
measurement basis is not the eigenstates of the density matrix, we have f n m{T, 0) ^ S n > m . 

In order to observe the effect of the measurement on the microscopic reversibility from a different viewpoint, let 
us derive the relation between pf(\ti(0)} — > \m(T))) and pn(Q\m(T)) — > 0|n(O))) in the operator-sum representation 
[l7j . The forward and the reversed transition probabilities are, instead of Eqs. Q and ((4]), expressed as 



MHO)) -> |m(T)» = Tr P m UP n p(0)P n WP m , (8) 
p R (9|m(T)) -> 0|n(O)» = Tr P n tf P m p(T)P m UP n , (9) 

n(0))(n(0)| and P m — \m(T))(m(T)\. Introducing the complementary operator Q n and Q m of P n and 



(k = n, m) 



> |m(T))) = TlP m UPnP(0)PnU f Pm, 

0|n(O)» = Tr P n U^P m p(T)P m Up n 

where P n 

Q k = i-P k , P k Q k = o : 

we can transform the forward transition probability as follows: 
M|n(0)) M T ))) = Trfr P m UP n p(0)P n 

= TrWP m p(T)(P m + g m )L>p- 1 (0) J PnP(0)P„ 
= Tr U^P m p(T)P m U(P n + Qn)M(0)P n p(0)P n 
+ Tr U^P m p{T)Qjjp- l {<S)P n p{<S)P n 

= Tr U^P m p(T)P m UP n p- 1 (0)PnP(0)Pn 

+ TrU^P m p{T)Qjjp- l {Q)P n p{Q)P n + TrU^P 



(10) 



Therefore, we have 



MKO)) ->• |m(T))) 
p R (9|m(T)) 9|n(0)» 



= 9n + Z, 



where 



9n 



(n(0)\p-\0)H0))(n(0)\p(0)H0)), 

Tr^tp mj5 ( r )Q m ^-i(Q)p w/5 (Q)p K + TrU''P m p{T)P m UQ n p- 1 (0)P n p(0)Pn 

p R (0|m(T)) 9|n(0))) 



(11) 

(12) 

(13) 
(14) 



We assumed that p(0) is the state where its inverse exists. In order to see the structure of the correction term £ more 
explicitly, let us divide it in the following way: 

£ = Pr 1 Trtf p m p(T)(Q m U + P m UQ n )M(0)P n p(0)Pn 
= Trtf P m p{T)[Q m U(P n + Q n ) + P m UQ n ]M{0)P n p{0)P n 
= Pr 1 Tr tfP m p(T)(Q m UP n + UQ n )p-\0)Pnp{0)Pn 
= Pr 1 g n Tr U^P mP (T)Q m Up n + p H 1 Tr P m p(T)UQ n p- l {Q)P n p(Q)P n 

+ (15) 
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As a special case, let us consider the situation where we measure the density matrix at t = with the eigenstate 
basis and obtain the result |n(0)), i.e. , p(0)P n = p n P n an d p~ 1 (0)P n = p~ x P n where p n is an eigenvalue of the density 
matrix p(0) for |n(0)). Then we have g n = PnPn 1 = 1 and £p — 0, but £ Q ^ £ 1Z: 

MHO)) \m(T))) 

+ (16) 



p R (e\m(T)) -> 9|n(0))) 
where 

C« = Pr 1 Tr U^P mP (T)Q m UP n . (17) 

Using ((5]), we can relate the correction term £ a to / nm (T, 0) as 

/„™(T,0) = 1 + £ Q . (18) 

When we measure the eigenstate of the density matrix at t = T, on the other hand, we have £ Q = 0. As we will 
explain in detail in Sec. lIVI for open systems, the corrections g n , £ a , and £^ depend on the choices of the measurement 
bases and, in the case of an open system, on the state of the total system as well. We can use g n ,£a and £ a as the 
measures of the irreversibility caused by the measurement processes. The analysis of (|12[) is more advantageous than 
flSJ) in the case of open quantum systems as we will see in Scc lIIII 



A. Example: One-dimensional free particle 



A free particle in one dimension is the simplest illustrative example. We assume that the initial state is given by 

P(0) = ||0i(O)><^(O)| + ||0 2 (O))<^ 2 (O)|, (19) 

where |0 fe (O)) (k 



1, 2) are the Gaussian wave packets in the position representation: 



(a#fc(0)) 



1 



(aH) 1 



(k) 

/4 ex P| l Po x 



X 

2^2 



(20) 



where x and p^ are the position and the momentum of the particle and the parameter a determines the width of 
the wave packet. Denoting M as the mass of the particle, we can write the state at t = T as 



p(T) = ±\MT))(MT)\ + l\MT))(MT)\, 



(x\MT)) 



iT 
Mm 



exp 



iT 
Ma 



2a 2 



ip' 'x 



2M 



(21) 
(22) 



We measure the initial state with the basis which contains |0i(O)) and 1 02 (0)) and measure the position x at t = T; 
we have g n = 1 and £p = 0. For the transition from 0i(O) to x{T), we have 

1 



1 



\(x(T)\MT))\ 2 + -\(x(T)\MT))\ 



(0! (0)1/3(0)10! (0)) 
(x(T)\p(T)\x(T)) 

MHO)) -+ |m(T)» - \(x(T)\U\M0))\ 2 (M0)\p(0)\M0)) 
which, according to (p~8|) . gives 

£ a = (\(x(T)\MT))\ 2 + \(x(T)\MT))\ 2 )- 1 



■MT)\MT))\ 2 , 



(23) 



(24) 



The T dependences of £ Q is plotted in Fig. [T] Here, we set p^ = 1, p {2> = 2, M = 1, and a = 1.2. The inverse of the 
correction term vanishes as x and T increase, which means that the correction grows much larger than unity; we 
cannot see the microscopic reversibility at all in this example. 

For the isolated quantum systems, we confirmed that the microscopic reversibility is not a general relation. Because 
there is nothing like thermalization, the deviation from the microscopic reversibility seems to depend sensitively on 
the choice of the system and the protocol; the microscopic reversibility is not a proper relation to characterize the 
dynamics of the isolated quantum systems. For the open quantum systems, nevertheless, we expect the microscopic 
reversibility is indeed a proper relation due to the effect of the thermalization; we will show in Sec. [V] that our 
expectation seems to be correct. 



FIG. 1: (Color online) The T dependence of £ a with = 1, p { 2) = 2, M = 1, and a = 1.2. 

III. MICROSCOPIC REVERSIBILITY IN OPEN QUANTUM SYSTEMS 

We next extend the discussion in Sec. |TT] to the case of open quantum systems. Let us consider the microscopic 
reversibility of the local quantum system which is thrown into a reservoir at t = 0. We allow the local system and the 
reservoir to be externally controlled by time-dependent parameters X s (t) and X r (t). The total Hamiltonian reads 

fltot(*) = H s (X s (t)) + £T r (Ar(t)) + H c 6(t), (25) 

where H s (X s (t)), H T (X T (t)), and H c are the Hamiltonians of the local system, the reservoir, and the coupling between 
the local system and the reservoir, respectively. The function 9(t) is a step function. Since we consider the situation 
where the coupling is turned on at t = 0, the initial state is given as a product state 

ptot(0) = Ps(0)<8p r (0). (26) 

As we mentioned in Sec. HI wc consider the forward and the reversed transition probabilities of the local system under 
the constraint that we measure the local system only. The transition probability of the forward process that the state 
of the local system evolves from |n(0)) to \m(T)) reads 

MKO)) -> \m(T))) = Tr r (m(T)|L>|n(0))(r i (0)|p t ot(0)|r i (0))(n(0)|L>t|m(T)) (27) 

and the reversed transition probability that the state of the local system evolves from 0\m(T)} to G)|n(0)) reads 

PR (Q\m(T)) -> 6|n(0)» = Tr r (n(0)|6>t|m(T))(m(T)|p tot (T)|m(T))(m(T)|L>|n(0)), (28) 

where (m(T)\U\n(0)} , (n(0)|p tot (0)|n(0)), and (n(0)|C7 t |m(T)) are the operators on the reservoir, whereas Tr r is the 
trace with respect to the degrees of freedom of the reservoir. As in the case of isolated systems, we again assumed 
that (m(T)\ep tot (T)Q\m(T)) = (m(T)\p tot {T)\m(T)) . Note that Pf(K°)) ~> \ m ( T ))) is indeed a real number since 
Ptot(0) is Hcrmitian: 

pp(|n(0)) ^ \m(T))y 

= ((rn(T), v\U\n(0), x) (n(0), x\p iot (0)\n(0), x') {n{Q),x'\U^\m{T),y) 

x,x',y 

= (m(T),y\U\n(0),x') (n(0), ^'1^(0)^(0), x) (n(0), x\fr\m(T), y) 

x,x' ,y 

= PF(|n(0)) -Hm(T)», (29) 
where x, x 1 , and y are the states of the reservoir. 
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As we did for isolated systems, let us introduce the projection operators P n to project on the state |n(0)) of the 
local system at time t — and Q n to project on the complementary space of P„, i.e. , 

P n = J2H0),x(0))(n(0),x(0)\, Q n = I-P n . (30) 

x(0) 

Similarly, we define P rn to project on the state \m(T)) of the local system at time t = T and Q m to project on the 
complementary space of P m , i.e. , 

P m = \m(T),x(T))(m(T),x(T)\, Q m = I - P m . (31) 

x(T) 

Now, we rewrite the forward probability as follows: 

Pf(KO)} -> \m(T))) = TrU^P m UP nPtot (0)P n 

= Tr ^tp m [/ i 5 tot (o)[7tc7 i 3 t -i( )p n ^ ot ( )p ra 

= Tr WP mPtot (T)Up^ t (0)P n p tot (0)P n 

= Tr U^P mPtot (T)(P m + Q m )J7/5 to 1 t (0)P„ / 5 tot (0)P„ 

= Tr L>t P mPtot (T)P m Up; o \ {0)P nPtot (0)P„ 

+ Tr L> t P m /5tot(P)Q m ^ to 1 t (0)P„/5tot(0)P„ 

(T)P m U(P n + Q n )p-\{0)P n p tot (0)P n 

+ Tr L> t P„ / 5 tot (T)Q m L>p t -i(0)Pn/5tot(0)P (1 
= Tr L> t P„ / 5 tot (r)P m L>P Il p t -t(0)P l /5tot(0)P l 

+ Tr U^P m ptot{T)P m UQ n p^(0)P n ptot(0)Pn 

+ TrU^P m p tot (T)Q m Upt o \(0)P n p tot (0)P n . (32) 



We assume that the initial state is the product state ()26|) . but the measurement bases are not necessarily the eigenstates 
of the density matrix of the local system, and thus the first term of Eq. (|32|) reads 



TrU^P m p tot (T)P m UP nP ^(0)P n p tot (0)P n = g nPR (e\m(T)) ->■ 9|n(0))), (33) 

where 

g n ee (n(0)\p;\0)\n(0))(n(0)\ Ps (0)H0)). (34) 

We do exactly the same transform as Eq. (fT5| for the second and third term except that the density matrix is /5 t ot 
and the projection operators are for the local system only. Then we arrive at one of the major results of the present 
paper: 

pp(|n(0)> ->■ |m(T))) 

- gn+Ka+Kfi, (35) 



p R (e\m(T)) -> 9|n(0))) 
where 

€a = Pr 1 9n Tr U*P m ptat {T)Q m U P n , 

& ee p^ 1 TrL>tp m p tot (T)^Q n(5to i( )p n/ 5 tot (0)p i . (36) 

We do not assume that the final state to be a product state. We cannot analyze the microscopic reversibility in the same 
way as in (fT2"j) in the case of open quantum systems, because the elements (n(0)|/5 to t(0)|n(0)) and (m(T)\p to t(T)\m(T)) 
are not c-numbcrs. Equation (|35|) reduces to (|12[) by eliminating the degrees of freedom of the reservoir. We will 
analyze the properties and the meanings of the corrections in Eq. (|35p in Sec. IIVI 
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A. A case where the initial and the final states are product states 

If the density matrix of the local system p s (0) in (|26|) is measured at t = and the result is an eigenstate \n(0)), 
we have p s (0)P n = p n P n , where p n is the eigenvalue of p s (0). Therefore, just as in the case of isolated systems, 
9n = Pn Pn = 1 an d £/9 = 0. If the final state is also a product state and the density matrix of the local system is 
measured at t = T with the result of an eigenstate |m(T)}, we have p s {T)P m = p m P m , where p m is the eigenvalue of 
p s (T); we have £ Q = in this case. When both of these conditions are satisfied, the microscopic reversibility holds 
exactly: 

pp(|n(0)) -> |m(T))) = m (Q\m(T)) -> 9|n(0))). (37) 

This is another major results of the present paper. This means that no matter how strong the local system is connected 
to the reservoir during the period between the measurements, the microscopic reversibility holds as long as the above 
conditions are satisfied. Note that this is a sufficient condition; we are not yet sure what is the necessary condition 
to make the correction terms vanish. 



B. A case where the local system is disconnected from the reservoir 

Again, we consider the case where we measure the density matrix of the local system at t = with the result of an 
eigenstate |n(0)). Now, let us consider the process where we gradually disconnect the local system from the reservoir 
and measure the energy of the local system at t = T with the result of an energy eigenstate |m(T)). Here we do not 
assume the form of the density matrix at t = T. 

If we split the time evolution at the time Ti so , at which we can regard that the local system is almost isolated from 
the reservoir, i.e. 

U(T) = U(T - r iso ) U(r iso ) = (cxp (-iH s (T - r iso )) <g> cxp (-iH T {T - r iso ))) U(n so ), (38) 

then we have 

pp(|n(0)) -> |m(T)» =Tr r (m(r)|L>(T iso )|n(0))(n(0)| / 5tot(0)|n(0))(n(0)|[/ t (T iso )|m(T)), 
PR (e\m(T)) -> 9|n(0)» = TrrMOjll^Ti^M^Xm^ (39) 

which arc independent of the state at t = T. Therefore, the transition probabilities and the correction term £ Q become 
constant after the local system is completely disconnected. 

Compared to the case of the isolated quantum systems, the behavior of the microscopic reversibility in open quantum 
systems is more nontrivial depending on the protocol that we choose. 

IV. DETAILS OF THE CORRECTIONS 

We only consider the case of open quantum systems because the consequences for isolated systems follow by 
eliminating the degrees of freedom of the reservoir. 

A. correction factor g„ 

Let us analyze the properties of the correction factor g n . Introducing a unitary matrix V that transforms the 
density matrix of the local system p s (0) in the present measurement basis to its diagonal form, i.e. 

Ps(0) = Vp D (0)V\ p D (0) = diag(pi,p 2 , • • ■ ,Pn), (40) 
we can write the general form of g n in (f3"4"]) as 

N 

9n = ~ \ Vni \ 2 l^'l 2 ' ( 41 ) 
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where Vij is the element of V. We have g n = 1 if we choose the state |n(0)) as an eigenstate of the density 

matrix of the local system. Since p s (0) is a positive operator, g n is positive, which is also obvious from (|41[) . The 
factor g n can be large if some states have a relatively small probability pj in the diagonalizing basis. Note, however, 
that g n is finite; to have g n = oo, it would require that pj = for the jth state, in which case the inverse of p s (0) 
would not exist. 



B. correction terms £ a and £p 

The correction terms from the microscopic reversibility must be quantities related to the disturbance due to the 
measurement process. Here, we consider the quantities a a = £, a pR and erg = £ppR. We recast these correction terms 
into the forms 

a a = g n TrP n U^P m p tot (T)(l - P m )U 

= g n Tr PJP P m p tot (T)U - g n Tr P n U^P m p tot (T)P m U 

= g n TvP n U^P m p tot {T)U -g nPR (e\m(T)) -> 9|n(0))), (42) 

and 

an = Tr P n Uip m p tot (T)U(l - P n )p t -i(0)P„Ptot(0) 

= TrP n frP m fa(T)U^(0)P n fa(0) -g n TrP n U^P m p tot (T)U 
= Tr P n U^P m UP n p tot (0) - g n Tr P n U^P m p tot (T)U 

= p F (|n(0)) -> |m(T)» - g n TrP n WP m p tot {T)U. (43) 
Taking the sum over the final states, we have from the third equality of (|4"2")) and the second equality of (|4^|) . 
a <* = HC^Tr^C/t^^^T);/ _ 9n Tr P n tf P m p tot (T)P m U) 

m rn 

= g n TvP n U^(J2 Pm)ptat(T)U - g n TrP n U^^2 p mPtot(T)P m )U 



= g n TrP n fr (p tot (T) P m p tot (T)P„}j U, 
X>/3 = ^(TrP n ?7 t P m ;7P„ i 5tot(0) - g n Tr P n W P m p tot (T)U) 

m m 

= Tr P n U^(^P m )UP n p tot {<S) -g n TrP n U^(J2P m )p tot (T)U 



(44) 



TrP„p tot (0) - 3„TrP n p to t(0) 

{l-g n )TrP n p tot (0). (45) 



The difference inside the parenthesis in (|44[) comes from the off-diagonal elements of the density matrix with respect 
to the degrees of freedom of the local system; therefore, the correction a a is the quantity related to the effect that the 
measurement at the final moment destroys the quantum coherence. The correction erg, on the other hand, depends 
on the factor g n and independent of the state at t = T; therefore, the quantity ap is related to the fact that the 
measured state at the initial moment differs from the eigenstate of the density matrix of the local system. The factor 
1 — g n indicates the difference from the eigenstate and TrP„/5 tot (0) is the probability of observing the state |n(0)). 



V. EXAMPLE OF AN OPEN QUANTUM SYSTEM: ONE-DIMENSIONAL SPIN CHAIN 



In order to evaluate the value of the correction from the microscopic reversibility quantitatively, we will numerically 
treat an open quantum system with a finite size reservoir. 
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A. Hamiltonian and protocol 

Let us consider the total system which consists of N pieces of 1/2-spins. We regard the first N s spins as the local 
system and the remaining N T {= N — N a ) spins as the reservoir. The Hamiltonian reads 

H tot (t)=H s (t)+H I (t)+H c 6[t), 

N B -1 

H s (t) = J2 J(S-S! +1 + 6(t)S?Sf +1 ) (46) 

i=l 
N a +N r -1 

H T {t)= Yl J(S?Sf +1 + 6(t)S?S? +1 ), 

i=N s +l 

H c = J (Sn s Sn b +i + Sn^a^+i) • (47) 

For t < 0, the total system consists of the two isolated Ising chains. We set them in the thermal equilibrium states 
at different inverse temperatures f3 s and /3 r . At t = 0, we measure the energy of the local system; the measurement 
basis of the initial state is the eigenstates of the density matrix of the local system, and thus g n = 1, <tq = erg = 0, 
and <jt = c Q . Then the local system is connected to the reservoir, so that for t > 0, the total system becomes 
an Xy-modcl. Finally, we again measure the energy of the local system at t = T. The strength of J is spatially 
uniform, and thus this is the case of strong coupling between the local system and the reservoir. This process satisfies 

the relation (m(T)| 9p(T)9|m(T)) = (m(T)\p(T)\m(T)) because the Hamiltonian is invariant under the flip of all 
spins. In order to calculate the time evolution, wc diagonalize the Hamiltonian of the total system numerically with 
LAPACK. 



B. Results 

We set N s = 2 and varied the number of the reservoir spins N T . Figured shows the T dependence of the forward 
transition probability pf(|«(0)) — > |?ti(T))) that we obtain the measurement results |«(0)) = | fj,) and \m(T)) = 
(| tt) + | ii))/V2, which are the energy eigenstates of the local system, while Fig. [2p shows the dependence of the 
corresponding reversed transition probability pR,(9|m(T)) — > 9|rt(0))). Figures and[2ji show the T dependences 
of the correction terms £ a and a a . We set the inverse temperatures j3 s = 1 and (3 r =0.1 with the coupling strength 
J = 0.1. As long as the number of the reservoir spins N T is finite, the finite-size effect appears as the measurement 
time T becomes large. We can, however, regard the local system as an open system up to T ~ 150 in the simulation 
for N = 10. It shows that the forward transition probability converges to a certain value which is presumably of a 
new equilibrium state for t > 0. The correction term £ Q converges to a very small (but finite) value, i.e., £ Q <C 1. It is 
rather surprising because, even though the local system approaches to an equilibrium state, the local system can be 
entangled with the reservoir strongly; the density matrix of the total system can be totally different from the product 
state and the measurement process could cause a large value of the correction term according to the discussion in 

secnni 

The correction seems to remain small generally in this model. Wc show in Fig. [3^, the same quantity as in Fig. [2}; 
in the case of N — 8, but by varying the value of the inverse temperatures of the local system (3 S and the reservoir 
f3 r . For T ~ 0, the ratio is sensitive to the value of the temperature of the local system /3 S . The ratio pf/pr for each 
parameter, however, seems to converge to a same small value as the local system goes to the new equilibrium state. 
Figure [3p shows the same quantity as in Fig. [2}; in the case of N = 8, but with the coupling strength J varied. The 
time evolution of the system is fast for the system with a large value of J, and hence the period for which the system 
indicates the behavior of the open system is short. In the region where we can regard the evolution of the ratio Pf/pr 
of each parameter as the behavior of the open quantum system (T < 90 for J = 0.08, T < 120 for J = 0.1, and 
T < 150 for J = 0.12), the ratios also seem to converge to a common small value as the local system goes to the new 
equilibrium state. Finally, Fig. 0^ and Fig. Hp show the cases of all the possible combinations of the states |n(0)) and 
|m(T)). Although the sign is negative in the case of |n(0)) = | Tt)j the absolute value of the ratio pf/pk of each case 
seems to converge to a common small value. It is difficult to determine from the numerical calculations whether all 
these small values of the ratio pf/pr at the stationary states coincide with each other, but they are of order 10~ 4 . 

Further theoretical study is required to estimate the order of the correction term compared to the transition 
probabilities. Nevertheless, the present simulation suggests that we can expect that the microscopic reversibility in 
open quantum systems almost holds even when the local system cannot be regarded as macroscopic. 
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FIG. 2: (Color online) The T dependence of quantities for the spin chain with ./V = 6, 7, 8, 9 and 10. The number 
of the system spins is iV s = 2. We set J = 0.1, /3 S = 1, j3 T — 0.1, and consider the transition from |n(0)} = ] ti) to 
= (I tt) + I 14)) A/2 as the forward process. The horizontal axis indicates the time T at which we measure the 
local system. The plots show the T dependence of (a) the forward transition probability pp(|n(0)) — > |m(T))), (b) the 
reversed transition probability pB.(Q\m(T)) — > O|n(0))), (c) the ratio of the forward and the reversed transition probabilities 
Pf/Pr = 1 + £<*, and (d) the correction term a a — £ a PR- 



VI. CONCLUSION 



We derived the correction terms of the microscopic reversibility of isolated quantum systems (0 and (fT2)) as well 
as of open quantum systems ()35j) by formal but exact treatment. Throughout the paper, we assumed the relation 

(m(T)| Q p(T)Q\m(T)) = (m(T)\p(T)\m(T)) and the product initial state for open quantum systems. 

We summarize the results of the present paper in TABLE I. For the microscopic reversibility in isolated quantum 
systems, we exemplified the case of a free particle system and found that the correction term can be very large. For 
the microscopic reversibility in open quantum systems, we first considered two situations which seem to be physically 
important: the case where the correction terms vanish (Sec. lIII A"|) and the case where we disconnect the local system 
from the reservoir in the middle of the time evolution (Sec. IIII B[) . In Sec. IIV1 we discussed the details of the 
corrections; we analyzed the bound of the factor g n and showed the meaning of the other correction terms explicitly 
by considering the quantities a a = J2 m C«Pr an d J2 m a P = Em m (03J an d (jUJ. Although we do not have 
an appropriate method of estimating the order of the correction terms theoretically, our numerical simulations of the 
one-dimensional spin chain suggested that, in the case of a thermal relaxation process, the correction term becomes 
very small compared to the transition probabilities; the microscopic reversibility almost holds even when the local 
system cannot be regarded as macroscopic. 

We expect that further analyses of the microscopic reversibility reveal more interesting properties of open quantum 
systems. 
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FIG. 3: (Color online) The T dependence of the ratio of the forward and the reversed transition probabilities in the case of 
N = 8 when (a) j3 s and /3 T are varied with J = 0.1 and (b) J is varied with /3 S = 1 and /3 T = 0.1. 




FIG. 4: (Color online) The T dependence of the ratio of the forward and the reversed transition probabilities in the case 
of iV = 8. (a) shows the result of all the transitions from |n(0)) = | ti) an d (b) shows the result of all the transitions from 
|n(0)) = | ft). As in Fig. H we set J = 0.1, & = 1, and /3 r = 0.1. 
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